Phase diagram for the Grover algorithm with static imperfections 
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We study effects of static inter-qubit interactions on the stability of the Grover quantum search 
algorithm. Our numerical and analytical results show existence of regular and chaotic phases depend- 
ing on the imperfection strength e. The critical border Ec between two phases drops polynomially 
with the number of qubits riq as Sc ~ nq '^^"^ . In the regular phase (e < Ec) the algorithm remains 
robust against imperfections showing the efficiency gain £c/e for e > 2~"'j/^. In the chaotic phase 
(e > Ec) the algorithm is completely destroyed. 

PACS numbers: 03.67.Lx, 24.10.Cn, 73.43.Nq 



the hardware properties of reahstic quantum computer 
0,0,0]. The effects of static imperfections on the ac- 
curacy of quantum computation have been investigated 
on the examples of quantum algorithms for the models 
of complex quantum dynamics [3, 0, Q, • 



suit a universal law for fidelity decay induced by static 



Quantum computations open new perspectives and 
possibilities for treatment of complex computational 
problems in a more efficient way with respect to algo- 
rithms based on the classical logic 0|. In the quantum 
computers classical bits are replaced by two-level quan- 
tum systems (qubits) and classical operations with bits 
are substituted by elementary unitary transformations 
(quantum gates). The elementary gates can be reduced 
to single qubit rotations and controlled two-qubit oper- 
ations, e.g. control-NOT gate Combinations of ele- 
mentary gates allow to implement any unitary operation 
on a quantum register, which for Uq qubits contains ex- 
ponentially many states N — 2"« . The two most famous 
quantum algorithms are the Shor algorithm for integer 
number factorization 01 and the Grover quantum search 
algorithm j3| • The Shor algorithm is exponentially faster 
than any known classical algorithm, while the Grover al- 
gorithm gives a quadratic speedup. 

In realistic quantum computations the elementary 
gates are never perfect and therefore it is very impor- 
tant to analyze the effects of imperfections and quan- 
tum errors on the algorithm accuracy. A usual model of 
quantum errors assumes that angles of unitary rotations 
fiuctuates randomly in time for any qubit in some small 
interval e near the exact angle values determined by the 
ideal algorithm. In this case a realistic quantum com- 
putation remains close to the ideal one up to a number 
of performed gates Ng ~ For example, the fidelity 

/ of computation, defined as a square of scalar product 
of quantum wavefunctions of ideal and perturbed algo- 
rithms, remains close to unity if a number of performed 
gates is smaller than Ng. This result has been estab- IV'(O) = C'lV'o) = sin ((f + l)a;G)|T)+cos ((i + 1)wg)|?7) j 

(1) 
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imperfections has been estabHshed for quantum al- 
gorithms simulating dynamics in the regime of quantum 
chaos. At the same time it has been realized that the ef- 
fects of static imperfections for dynamics in an integrable 
regime are not universal and more complicated. There- 
fore it is important to investigate the effects of static 
imperfections on an example of the well known Grover 
algorithm. First attempt was done recently in 0|, but 
the global picture of the phenomenon remained unclear. 
In this paper we present extensive numerical and analyt- 
ical studies which establish the global stability diagram 
of reliable operability of the Grover algorithm. 

Let us first outline the key features of the Grover al- 
gorithm 0. An unstructured database is presented by 
N = 2"9 states of quantum register with Uq qubits: {|a;)}, 
a:: = 0,...,A^— 1. The searched state |t) can be identified 
by oracle function g{x), defined as g(x) = 1 if x = r and 
g{x) = otherwise. The Grover iteration operator G is 
a product of two operators: G — DO. Here the oracle 
operator O = (— 1)9(^) is specific to the searched state 
|r), while the diffusion operator D is independent of |r): 
Du = — 1 -I- and D^j — ^ (* 7^ j)- For the initial 
state 1-00) = J2x=^ t applications of the Grover 

operator G give [l|: 



lished analytically and numericallv in extensive studies 
of various quantum algorithms jj, lE IE S S IE 0|- 

Another source of quantum errors comes from inter- 
nal imperfections generated by residual static couplings 
between qubits and one-qubit energy level shifts which 
fiuctuate from one qubit to another but remain static in 
time. These static imperfections may lead to appearance 
of many-body quantum chaos, which modifies strongly 



where the Grover frequency ujg — 2 arcsin(-y/l/A^) and 

\ri) — X^i'^T^^^'' \^) /VN — 1. Hence, the ideal algorithm 
gives a rotation in the 2D plane (|r), |?7)). 

The implementation of the operator D through the el- 
ementary gates requires an ancilla qubit. As a result the 
Hilbert space becomes a sum of two subspaces {\x)} and 
{|a; -I- iV)}, which differ only by a value of {ug + l)-th 
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qubit. These subspaces are invariant with respect to op- 
erators O and L»: O = 1 - 2\t){t\ - 2|t + N){t + N\, 
D = l- 2|^o>(V'o| - 2|^i)(^i|, where = 1^ + 

N)/^/N and iV'o.i) correspond to up/down ancilla states. 
Then D can be represented as D = WRW where 
the transformation W = Wn^ ■ ■ ■ Wk ■ ■ ■ W\ is composed 
from Tiq one-qubit Hadamard gates , and R is the n^- 
controlled phase shift defined as Rij —Oiiiy^j, Rqo = 1 
and Rii = — 1 if i 7^ (z, j = 0, . . . , iV — 1). In turn, this 
operator can be represented as i? = Wn^f^n^-i ■ • ■ ""i 
(J^^_i . . .afWrig, where is generalized n^-qubit Tof- 
folli gate, which inverts the rtg-th qubit if the first Uq — l 
qubits are in the state |1). The construction of An, from 
3-qubit Toffolli gates with the help of only one auxil- 
iary qubit is described in [17]. As a result the Grover 
operator G is implemented through Ug = 12ntot — 42 
elementary gates including one-qubit rotations, control- 
NOT and Toffolli gates. Here ritot = -I- 1 is the total 
number of qubits. 

To study effects of static imperfections on the Grover 
algorithm we use the model introduced in In this 

model a quantum computer hardware is described by the 
Hamiltonian H: 

i i i<j 

Here, ct^ are the Pauli matrices for qubits i, and A is 
an average one-qubit energy spacing. All ritot qubits are 
placed on a rectangular lattice, the second sum in Hs 
runs over nearest neighbor qubits with periodic bound- 
ary conditions. Qubit energy shifts Ui and couplings bij 
are randomly and uniformly distributed in the intervals 
[—a, a] and [—13, f3], respectively. Following 
we assume that the average spacing A is compensated by 
specially applied laser pulses so that between subsequent 
elementary gates the wavefunction evolution is given by 
the propagator Us = exp{—iHstg). Thus all static er- 
rors are expressed via this propagator while the elemen- 
tary gates are taken to be perfect. Appropriate rescaling 
of parameters and 5^ allows to put tg = 1 without 
any loss of generality. We concentrate our studies on the 
case a — (3 = e where inter-qubit couplings lead to a 
developed quantum chaos |llLll4l|. 

A typical example of imperfection effects on the accu- 
racy of the Grover algorithm is shown in FigQlfor a fixed 
disorder realization of Hs in (|2ll on 3 x 4 qubit lattice. 
It clearly shows that imperfections suppress the proba- 
bility WG to find the searched state, where wg is given 
by a sum of probabilities of states |t) and |t -|- N). In 
contrast to the case of time-dependent random quantum 
errors studied in in the case of static imperfections the 
oscillations of probability wg do not decrease with time 
t. Another interesting feature is a significant decrease 
of the period of the Grover oscillations compared to the 
case of ideal algorithm, where Tg = 7r/2tJG. This effect 
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Figure 1: Probability of searched state wait) (top) and fi- 
delity f{t) (bottom) as a function of the iteration step t in 
the Grover algorithm for ritot = 12 qubits. Dotted curves 
show results for the ideal algorithm (e = 0) , dashed and solid 
curves correspond to imperfection strength e — 4 ■ 10~^ and 
10~^, respectively. 

is also absent in the case of random errors. The fidelity of 
quantum computation f{t) also shows non-decaying os- 
cillations at large times. However, in average the maxima 
of fidelity correspond to minima rather than maxima of 
probability wg- Hence, f{t) is not appropriate for tests 
of the algorithm accuracy. 

Following a pictorial presentation of the dynam- 
ical evolution in the Grover algorithm can be obtained 
with the help of the Husimi function which is shown 
in FiglSl In this presentation the computational basis 
X can be considered as a coordinate space representa- 
tion for the wavefunction %px{t) {x = 0,...,2N — 1), 
while the conjugated basis obtained by the Fourier trans- 
form corresponds to momentum representation p {p — 
— iV-l-1, . . . ,N). In this way the initial state of the Grover 
algorithm \'ipo) gives a peaked distribution with p — 0. In 
the ideal algorithm the total probability is distributed be- 
tween two states |r) and \r]) (see Eq.l^) that gives two 
orthogonal lines in the phase space of Husimi function 
(see FiglSl top raw). After the period Tg ~ 34 all the 
probability is transferred to the target state |r) {wg ~ 1). 
In the presence of moderate imperfections the fiips of the 
ancilla qubit become possible that involves into dynam- 
ics two additional states. As a result the probability is 
mainly distributed over four states corresponding to four 
straight lines in phase space (Figl^l middle raw): 

ko>-|r) \n) = \T + N) 

The probability W4 contained in these states is close to 
unity (in Fig|2| W4 = 0.998 for e = 10^^). Above cer- 
tain critical border Sc this simple structure is completely 
washed out {w4 — 6 • 10^''), and the Husimi function 
shows only random distribution (Fig[21 bottom raw). 
The dominant contribution of these four states can be 
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Figure 2: (Color) Evolution of the Husimi function in the 
Grover algorithm at times t — 0, 17, and 34 (from left to 
right), and for e = 0, 0.001, and 0.008 (from top to bottom). 
The qubit lattice and disorder realization are the same as 
in FigH The vertical axis shows the computational basis 
X = 0, . . . , 2N — 1, while the horizontal axis represents the 
conjugated momentum basis. Density is proportional to color 
changing from maximum (red) to zero (blue). 



also seen in spectral density S{uj) of the wavefunction 
ipx{t)- This density is defined as: S{Ld) — |a^(w)p, 
where axiuS) = ^^Zq exp(iu;i)/-y/T/ and T/ is a 
large time scale on which the spectrum is determined 
(we usually used T/ ~ 5Tg ^ Tg). The phase diagram 
of spectral density S{uj) dependence on the imperfection 
strength e is shown in FigEl Two phases are clearly seen: 
for e < Ec the diagram contains four lines corresponding 
to the four states while for e > Sc these Hnes are 
destroyed and the spectrum becomes continuous. These 
phases correspond to the qualitative change of the Husimi 
distribution shown in FigEl 

To study the transition between these phases in a more 
quantitative way we analyze the dependence of probabil- 
ities wq and W4 on the imperfection strength e for a large 
number of disorder realizations in Hs lO changing also 
the number of qubits ritot- The number of reaHzations 
vary from 50 to 1000 depending on e and ritot- Since 
the frequency of Grover oscillations varies strongly with 
e and disorder we average wg and over a large time in- 
terval Tf to suppress fluctuations in time. The obtained 
results are summarized in Fig^ For a fixed value of ntot 
the dependence wg{£) changes strongly from one reaHza- 
tion to another (Fig^). In contrast, the probability wa 
remains close to unity being insensitive to variations of 
disorder up to £ < £c (Fig^b)- Only for e > £c, when 
W4 ^ 1, it becomes sensitive to disorder. The proba- 
bilities averaged over disorder wg and W4 are shown in 
FigBi,b. They also have a qualitative change in behav- 
ior near £c, especially W4. These results confirm the fact 
that the phase transition takes place near some critical 
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Figure 3: (Color) Phase diagram for the spectral density S{uj) 
as a function of imperfection strength e, ntot = 12, same dis- 
order realization as in FigEl Color is proportional to density 
S{ui) (yellow for maximum and blue for zero). 



£c for an ensemble of disorder realizations. 

The value of £c can be obtained from the following 
estimate. The transition rate induced by imperfections 
after one Grover iteration is given by the Fermi golden 
rule: F ~ e'^n'^ntot, where ntot appears due to random 
contribution of qubit coupHngs e while factor takes 
into account coherent accumulation of perturbation on 
Ug gates used in one iteration (see, e.g. In the 

Grover algorithm the four states ^ are separated from 
all other states by energy gap AE ~ 1 (it appears due 
to sign change introduced by operators O and D). Thus 
these four states become mixed with all others for 



£ > £c 



(4) 



when r > AE. Here the numerical factor is obtained 
from numerical data. The parameter dependence is well 
confirmed by data for u>4 shown in Fig^li. 

The variation of averaged Grover probability tmg with 
£ and ntot is presented in Fig^t- The dependence on 
system parameters can be understood on the basis of 
simple single-kick model. In this model the action of 
static imperfections in all gates entering in one Grover 
iteration is replaced by a single kick unitary opera- 
tor Ueff = exp {—iHsngR) acting after each iteration. 
Here i? is a dimensionless renormalization factor which 
takes into account that gates do not commute with Hs. 
Figs3K)b show that this single kick approximation gives a 
good description of original averaged data with R — 0.56. 
Thus, the renormaHzation effects play a significant role 
and therefore this model does not describe the probabil- 
ity variation for a given disorder realization. However, 
the averaged dependence is correctly reproduced. 
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Figure 4: (Color online) Dependence of probabilities wa (a,c) 
and Wi (b,d) on rescaled imperfection strength e/sc, with £c 
from Q. For panels (a,b) ritot ~ 12, squares and pluses 
show data for two typical disorder realizations, green/grey 
area shows the region of probability variation for various dis- 
order realizations (see text), full thick curves give average de- 
pendence wg, W4. Dashed area bounded by thin curves show 
the region of probability variation in the single-kick model, 
open circles give the average data in this model with rescal- 
ing factor R = 0.56. Panels (c,d) show wc, W4 for ntot = 9 
(triangles), 12 (full circles), 15 (open squares) and 16 (full 
squares). In panel (c) full curves are given by Eq.J^ for same 
ntot values from top to bottom, R = 0.56. 



In the regime where the dynamics of Grover algorithm 
is dominated by four states subspace lO the single-kick 
model can be treated analytically. The matrix elements 
of the effective Hamiltonian in this space are 



H, 



/ A + a 




eff 



lUJG 





-iujG ^ 

A — a —ii^G 

B b 

ilog b B 



(5) 



where A = —Rug o-ii'^'W^ 



(^)i 



\T),B = RngJ2:<j-i,j 
b, a = -Rrigan^+i and b = Rng{bn^+i,n,+2-L^ + 
&n,+i,L, + &n,,n,+i + &n,+i-L, ,«,+! ) and QuMts are ar- 
ranged on Lx X Ly lattice, and numerated &s i = x + 
Lx{y - 1), with X = l,...,Lx, y = l,...,Ly. In the 
limit of large Uq the terms a, b are small compared to 
v4, i? by a factor 1 / y/n^ and H^f / is reduced to 2 x 2 
matrix, which gives wg = 2uj'q/[{A — B)"^ + 4cJq]. For 
large Uq the difference A — B has a Gaussian distribution 
with width a — RngyJnq/?>^J -\- 2/3^ = eRug^Jn^. The 
convolution of wg with this distribution gives 

WG = \/t^{1 -erf(\/2tJG/cr))exp (2^^/0-2) ujg/ct (6) 

This formula gives a good description of numerical data 
in Fig^t that confirms the validity of single-kick model. 
For a S> WG and a typical disorder reaHzation with 



{A — B) ^ a the actual frequency of Grover oscillations 
is strongly renormalized: w « (A — _B) ~ cr ^ log, and 
in agreement with FigElw ~ e/Sc- In this typical case 
Wg ~ ^'gI<^'^ ^ 1 (almost total probability is in the 
states |'7o),|'7i)). Hence, the total number of quantum 
operations Nop, required for detection of searched state 
|r), can be estimated as Nop ~ Nm/^ ^ '^I^'g ^ ^^/^c, 
where Nm ~ ^ '^^/^g is a number of measure- 

ments required for detection of searched state jl^l. Thus, 
in presence of strong static imperfections the parametric 
efficiency gain of the Grover algorithm compared to clas- 
sical one is of the order £c/e. For e ~ log the efficiency is 
comparable with that of the ideal Grover algorithm while 
for e ^ Ec there is no gain compared to the classical case. 

In summary, we have shown that the Grover algorithm 
remains robust against static imperfections inside a well 
defined domain and determined the dependence of algo- 
rithm efficiency on the imperfection strength. 
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